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^ ■ There have been suggestions to measure atomic parity nonconservation (PNC) 

atomic structure effects. Precise atomic PNC measurements could make a significant 
Q> ■ contribution to tests of the Standard Model at the level of one loop radiative correc- 



tions. However, the results also depend upon certain features of nuclear structure, 
such as the spatial distribution of neutrons in the nucleus. To examine the sensitivity 
to nuclear structure, we consider the case of Pb isotopes using various recent rela- 



^ ■ tivistic and non-relativistic nuclear model calculations. Contributions from nucleon 

internal weak structure are included, but found to be fairly negligible. The spread 
among present models in predicted sizes of nuclear structure effects may preclude 
using Pb isotope ratios to test the Standard Model at better than a one percent 
level, unless there are adequate independent tests of the nuclear models by various 
alternative strong and electroweak nuclear probes. On the other hand, sufficiently ac- 
curate atomic PNC experiments would provide a unique method to measure neutron 
distributions in heavy nuclei. 
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I. INTRODUCTION 



Parity nonconservation (PNC) in atoms arises from the electroweak interaction between 
the electrons and nucleons, primarily due to exchange of the neutral gauge boson, Zq. The 
dominant contribution in heavy atoms comes from the coupling of the axial electronic current 
to the vector nucleon current. Because the vector currents are conserved, atomic PNC 
essentially measures the electroweak coupling to the elementary quarks, bypassing many 
of the difficulties of hadronic physics. Thus in principle atomic experiments can measure 
certain key electroweak parameters quite accurately, and also help probe for new physics 
beyond the presently successful Standard Model of the electroweak interactions. 

In fact, there remains much to be learned about the Standard Model, including the 
masses of the top quark and the predicted Higgs boson(s), and whether there are additional 
generations of quarks and leptons. In addition the Standard Model faces the well-known 
gauge hierarchy problem, and it is certainly possible that electroweak measurements may 
reveal something totally new, such as technicolor or supersymmetric particles. Accurate 
measurements of PNC in atomic cesium already play an important role in addressing such 
questions. 

Two major issues affect the interpretation of atomic experiments and will become more 
crucial as experimental accuracy improves, namely the small but not negligible effects of 
nuclear size and structure, ^ and the reliability of the atomic theory of heavy atoms 
Atomic theory, the source of the largest uncertainty, has received a great deal of attention 
leading to increasingly precise calculations of PNC for a number of elements. Cesium, in 
particular, is now believed to be understood at the 1% level. 

To advance further may require canceling out all uncertainties of the atomic theory by 
comparing PNC measurements on different isotopes of the same element. Such experiments 
in fact have been proposed using strings of isotopes of such elements as Cs, Dy, and Pb. 

As we discussed in a previous note, hereafter referred to as I, it then becomes im- 
portant to find the level at which nuclear structure interferes with interpreting atomic PNC 
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purely in terms of particle theory. The wave function of the atomic electrons varies over 
the dimensions of the nucleus, causing the net electroweak interaction with the nucleons to 
depend on the spatial distribution of both the protons and neutrons. As demonstrated in I, 
the PNC observable is (for sin^ 9w ~ 1/4 and Rn ~ Rp) very roughly proportional to 

l-A(Zani + 5</i?g + ..., (1.1) 

where Rn and Rp are the equivalent rms radii for the nuclear distribution of neutrons and 
protons. The proton, or rather the nuclear charge, distribution is well known from elec- 
tric probes: electron and muon scattering, optical isotope shifts, muonic atoms, etc. The 
extraction of the neutron distribution, however, is quite model-dependent and difficult to 
determine to the same high accuracy. 

On the one hand, the neutron distribution is needed in order to extract the weak param- 
eters in heavy atom experiments. On the other hand, to the extent the weak parameters 
are known, the experiments provide a method of measuring the changes in the neutron dis- 
tribution, primarily the rms radius. Thus atomic experiments on isotopes of heavy atoms 
may provide a unique opportunity to test nuclear model calculations. 

In this paper we explore the nuclear structure issues extensively. We have utilized several 
recent detailed nuclear structure calculations from various authors, in order to quantitatively 
estimate the nuclear model-dependent corrections to atomic PNC in the isotopes of 82Pb, 
an element of interest experimentally. Our major conclusions are: 

(1) For single isotope measurements on ^''^Pb, (Za) is sufficiently large that the effects 
of nuclear structure on atomic PNC cannot be neglected. The uncertainties due to neutron 
distributions appear to be less significant for the extraction of electroweak parameters than 
the current uncertainties due to atomic structure. However, unless the change in neutron 
distributions along an isotopic chain can be better predicted (or independently measured, 
e.g. via parity violating e~ scattering), PNC ratios of 82 Pb isotopes will not be able to 
provide an extraction of the weak mixing angle to better than a one percent level. 

(2) In lighter nuclei (including the important case of 55 Cs), Za is sufficiently small 
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that uncertainties in nuclear structure can probably be safely ignored compared to current 
uncertainties in atomic theory, when using a single isotope. Calculations for the non-magic, 
odd-Z, Cs nuclei pose additional difficulties, and further investigation is still required to 
determine how accurately the Standard Model could be tested when using ratios of isotopes. 

(3) Atomic PNC experiments provide perhaps one of the cleanest opportunities yet avail- 
able to study the nuclear neutron skin. The situation is similar to atomic isotope shifts which 
have provided precise measurements of the ratios of changes in the nuclear rms charge radius 
in strings of isotopes. Here we have a weak probe of the neutron distribution which is free of 
the gross uncertainties associated with strongly interacting probes. At the level of precision 
that the Standard Model is known, this yields another testing ground for nuclear models. 

We note that alternative weak probes, such as parity-violating intermediate-energy e"- 
nucleus cross section asymmetries, are also sensitive to the neutron distributions. In combi- 
nation with atomic PNC, these may help simplify the separation of nuclear structure effects 
from electroweak radiative corrections. 

The paper is organized as follows: In section II we sketch briefly the relevant parts 
of electroweak theory. In section III we review the simple analytical model we presented 
in reference (1); this provides a convenient framework for discussing the effects of nuclear 
structure on atomic PNC in terms of the rms radii of the proton and neutron distribu- 
tions in the nucleus. In section IV we discuss the electroweak nucleon form factors. The 
intrinsic electroweak structure does begin to contribute at the level we are interested in, 
but the uncertainties in this structure (due to strangeness admixtures, etc.) should have 
a negligible effect on the total PNC amplitude of the nucleus. In section V, we examine 
the key ingredients of currently available theoretical models for heavy nuclei, including both 
non-relativistic and relativistic formalisms. We consider the reliability of these models, and 
discuss the need for calculations which include correlations beyond the Hartree-Fock level. 
We also discuss alternative experimental means to measure the desired neutron distribu- 
tions. In section VI we discuss the relevant Standard Model parameters, and the accuracy 
desired in their extraction from atomic parity violation. We derive the propagation of error 
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from nuclear model uncertainties to electroweak parameters, focusing on isotopes of lead 
and (to a lesser extent) cesium. These two elements are of current experimental interest, 
and are representative of the very heavy and moderately heavy regions of the periodic table. 
In section VII we discuss our numerical results, using various existing Hartree-Fock nuclear 
calculations, summarize, and discuss the conclusions from the previous sections. 



Because the Z-boson is massive (91.16±.03 GeV), the quark-electron interaction due to 
Z-exchange may be taken to be of zero range compared to atomic or nuclear dimensions. 
What we observe in atoms is the electron interaction with nucleons, not individual quarks. 
Nucleons are, of course, composite structures consisting each of three quarks net, but also qq 
pairs as well as gluons. We make the assumption here, which we justify in Sec. IV, that we 
can neglect the internal nucleon structure and simply add the point coupling of the 3 quarks 
to obtain the net nucleon weak couphng. The PNC part of the nucleon-electron interaction 
can be written in terms of axial and vector currents 



If, in addition to neglecting internal nucleon structure, we treat the nucleons nonrelativisti- 
cally (a very good approximation), we have 



II. PNC IN THE STANDARD MODEL 



A. Theoretical Considerations 



HpNC = Vn X Ae + An xVe. 



(2.1) 





(2.2a) 



where B stands for n (neutron) or p (proton) and 
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C2p = l{l-4:sm'^ew)gA, 

Cm = -|, (2.2b) 
Can = -1(1 -4:8110^ 9w)g A- 

These expressions assume tree-level Standard Model couplings. In section VI, we discuss 
the important effects of radiative loop corrections. 

The first term in eqn. (2.2a) grows coherently with nucleon numbers and Z. The 
second term, together with the anapole moment [^] term (which also depends upon ctb ■ a), 
amounts to at most a few percent of the first term in heavy atoms, and furthermore sums 
to zero when all hfs sublevels are combined, since all directions of (Tb are then weighted 
equally. Thus in this paper we will consider the first term only. The effective interaction is 



HpNCi = ^J [-Npn{r) + Z{1-A sin^ ew)ppir, 

X il)l'^^%l)^d\ , (2.3) 



where here the p„ and pp are normalized to unity. The neutron and proton densities include 
a folding with the weak form factors (see Sec. IV). 

We need the spatial variation of ipl'-f^ipe over the nucleus, its normalization, and its 
dependence on nuclear structure. PNC effects are dominated by Si/2-electrons (k = —1) 
coupled to pi/2-electrons (k = +1). We define 

p,{r) = tPlin-f'tsin , (2.4) 

which turns out to depend only on the magnitude of r. P5(r) can be factored conveniently 
as follows: 

p,{r) = C{Z)Af{Z,R)f{r), (2.5) 

where C{Z) contains all atomic structure effects for a point nucleus including many-body 
correlations; Af = il'l{^)l^'^s{^) is the normalization factor for a single electron; /(r) contains 
the spatial variation and is normalized to /(O) = 1. 



Because the electric potential is very strong near the nucleus, we can safely neglect atomic 
binding energies in /(r). In Pb, for example, the potential at the nuclear surface is about 
15 MeV compared with valence electron binding energies of a few eV. In addition, to a very 
good approximation, |jT[] 

M = R-\ (2.6) 



where 7 = 2 1 — yl — {ZaY and i?, often called the equivalent charge radius, is given by 

nl/2 



R 



2 < r ^charge 



(2.7) 



We are not interested in the absolute value of A/", but only its variation with nuclear structure. 
Observable PNC effects are proportional to the matrix element between two atomic states 
i and j 



G 



<i\Hpj,c,i\3>= ^C.,,{Z)N 



2V2 



X 



[- Nqn + Z{l-Asm' ew)qp]. 



(2.8) 



As mentioned above, this is modified by radiative corrections which we discuss in some detail 
in section VI. Effects of nuclear structure on PNC are contained in M and the two quantities 



qn = J Pn{r)f{r)d'r, (2.9 a) 

Qp = f Pp{r)f{r)d\. (2.9 b) 



We note that 1 — 4 sin^ 6^ is a small number; from high energy experiments, sin^ 6w = 
0.230 ± .004. The value of sin^ 9w can also be deduced from atomic experiments with an 
accuracy that will be limited in part by nuclear structure effects, as we discuss in sections 
VI and VII. 

The proton (charge) nuclear form factors needed for and JV are generally well known 
from measurements of the charge distribution of nuclei close to the stable valley and many 
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unstable nuclides as well. Neutron nuclear form factors are needed for g„, and are not well- 
determined experimentally, and statements about them are quite model-dependent. Neutron 
and proton distributions are often taken to be proportional to each other, scaled by and 
Z. However, neutron-rich nuclei have larger neutron distributions than the protons and the 
reverse is true for proton rich nuclei. In an isotopic sequence, the A^l'^ law is not followed 
for either the charge or the neutron distributions separately. 



B. Experimental Considerations 

As first pointed out by Bouchiat and Bouchiat, 0] the effect of Hp^c in neutral atoms 
grows rapidly with atomic number Z, approximately as Z^. Thus experimental interest has 
concentrated on heavy atoms, namely 55CS, siTl, 82Pb, and ssBi. (For some reviews, see 
reference |^.) The measured quantity in all experiments is the electric dipole amplitude 
£pNC between two electronic states which, in the absence of HpNc, would have the same 
parity and hence would have no electric dipole amplitude connecting them. Denoting the 
initial and final states by i and /, we can write: 

■< /|Ei|r7, >< nlHpNcli > 



£pNc — 



+ 



< f\HpMc\n X n\Ei\i > 



(2.10) 



Wf-Wn 

where the first and second terms give the mixing due to Hp^c of opposite parity states 
into the initial and final states respectively. W is the energy of the atomic states, and 
El = — I]jerj is the electric dipole charge operator. The magnitude of £pNC is of order 
lO^^eao for the heaviest atoms of interest. 

Two experimental techniques have evolved for measuring £pnc- One involves applying 
an external static electric field which, like Hpjyc, mixes in opposite parity states and creates 
an electric dipole amplitude between the states i and /. The interference between this 
Stark amplitude and SpNC leads to a parity- violating signature in the optical transition 
from i to / in which the sign of the interference term reverses with the sense of circular 
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polarization of the incident light, and with other vectors specifying the handedness of the 
experimental arrangement. The other technique uses no external fields, but instead exploits 
the interference between £pNC and the magnetic dipole (Ml) amplitude between the same 
two states. This interference causes parity-violating optical rotation, i.e. a rotation of the 
plane of polarized light passing through the atomic vapor at wavelengths near the magnetic 
dipole absorption line. 

The Stark interference technique has been used in the all-important measurement of 
PNC in Cs on the highly forbidden QSi - 7S^ Ml absorption line at 532 nm, and in the 
measurement on the 6Pi - 7Pi Ml transition in Tl. The optical rotation technique has been 

2 2 

applied to the allowed Ml absorption lines at 876 nm and 648 nm in Bi, and to the similar 
1278 nm and 1283 nm lines in Tl and Pb respectively; all of which involve transitions among 
low-lying configurations of 6p electrons, for example 6Pi - 6P3 in Tl. Both techniques have 
reached the 1 percent level of accuracy. 

Among the elements studied thus far, Cs and Pb are the most hkely candidates for 
comparing different isotopes. It may be possible in the case of Cs to use optical atom traps 
to carry out measurements on a long string of radioactive isotopes. Measurements on Pb 
will probably be restricted to stable or long-hved isotopes. In either case, achieving the 
level of accuracy discussed in this paper (a few percent down to 1 percent in the isotopic 
difference) , although possible in principle, will be a challenging task in the next generation 
of atomic PNC experiments. 

III. A SIMPLE MODEL FOR THE NUCLEAR FORM FACTORS 

Given proton and neutron distribution functions, there is no difficulty in calculating 
gp, Qn and the variation in M. In I, we used a simple model to estimate the importance of 
nuclear structure on PNC observables. We review those results here. 

Consider a uniform nuclear charge distribution of radius R. This charge produces an 
electric potential 
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^ ' l-l/r, r> R. ^ ' 

A power series for the Dirac wave function inside the nucleus yields 

/(r) = 1 - i(Za)^[(r/i?)^ - \{r/RY + ^(r/i?)^ 

+ (3.2) 

Again for the sake of simplicity here, we assume that, as for a uniform distribution (for 
either n or p), < >= and < >= where i?^ = | < >. From (2.9a) we find 
(neglecting here any differences between charge and proton radii) 

qp^l-Q.2m{Zaf + 0{Zaf , (3.3 a) 

which is insensitive to nuclear structure to this order. From (2.9b), we find 

^ 1 - {Zaf 1^0.038 + 0.221^j + O^Zaf , (3.3b) 

which does depend on the neutron form factor. Here we have introduced equivalent neutron 
and proton radii of the form (2.7); the second form in (3.3b) assumes that (Rn/Rp)^ — 1 is 
small. 

In this section, we have made rough approximations in order to illustrate the sensitivity 
of the results to moments of the neutron and proton distributions. For comparison with 
experiment, a more detailed analysis is necessary, using actual solutions of the Dirac equation 
for realistic charge distributions and the best available theoretical neutron distributions. 
This is done in Sees. VI and VII. 

IV. INTRINSIC NUCLEON STRUCTURE EFFECTS 

The usual treatment of atomic PNC begins with an effective Hamiltonian for the par- 
ity violating electron-nucleus interaction, as in equation (2.2a), which involves normalized 
proton and neutron distributions 
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Zpp{r) = (4.1a) 
p 

Npnin = T^i^anMr)), (4.1b) 

n 

where ijj^jp is a destruction (creation) operator for nucleons, and the matrix elements are be- 
tween nuclear ground states. However, these formulae implicitly assume point-like nature for 
nucleons, and thus the usual analysis makes no distinction between weak, electromagnetic, 
or point nucleon distributions, aside from overall charges. 

Of course, nucleons do have an internal structure, and this must be properly folded into 
the above distributions. The internal weak structure is related to, but different from, the 
electromagnetic structure, and can be calculated in the context of the Standard Model. We 
demonstrate in this section that the known, electromagnetic structure of nucleons yields a 
rather small overall effect on atomic PNC calculations, but must be included when extremely 
high precision results are required. 

There has been considerable discussion in recent literature P,|9|,p!0| concerning the pos- 



sibility of nontrivial strange quark matrix elements in the nucleon. This could lead to a 
sizable "strangeness radius" of the nucleon, which in turn would modify the weak radius in 
a well defined way. We allow for this possibility in our analysis, although such a strangeness 
contribution to atomic PNC is likely to be quite negligible. 

In the Standard Model, assuming in addition that strong SU(2) isospin is a good sym- 
metry for the nucleons, one can extract relations between weak and electromagnetic form 
factors inil which then describe the internal nucleon structure: 



GT'^^iq') = '^{1- Asm' ew)GY{q') 

-l{Gr{q') + GUq')), (4.2a) 
= i(l-4sin2^v^)Gr(g') 

-k{Griq') + GUq'))- (4.2b) 

Here, G^'^ is the usual Sachs electric form factor for a current operator , where X can 
represent weak, electromagnetic, or specific quark flavor currents: 
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<p\N\J^\p,N> 

+zF2-^'^a^,g7(2M))M(p), (4.3a) 
GW) = F^{q') + {qV^M')F^{q'), (4.3b) 



and 

= |(m7^m) - i(J7^rf + S7^s) (4.4a) 

= (l-fsin^MTO 

+ (-| + ism^M(^,' + ^^), (4.4c) 

are the Standard Model electromagnetic and weak vector currents in terms of quark field 
operators. (We ignore quarks heavier than strange.) C^^ is thus the strangeness electric 
form factor, and is constrained to be strictly at = 0. Note that one recent estimate 
gives a strangeness mean-square-radius of around 0.14 fm^, roughly as large as that for the 
neutron electric charge (but of opposite sign). (This quantity can in principle be measured 
in, e.g. parity violating scattering from nucleons at forward angles.) 

With the above relations, we see immediately that at = 0, the usual weak charges are 
exactly obtained: 

Q; = i(l-4sin2^vF), (4.5a) 
Q: = -h (4.5b) 
and one can also predict weak rms radii 



2 \' // 



+ 6(i)(l-4sin2 ew)/{8M^), 



(4.6a) 



I,n 2^ J \ I l^n 2 { ^/jP 



rT = ia-4sin2^^)(rn" - i 



-|(^'). + 6(-i)/(8M2). (4.6b) 
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where the subscript I indicates intrinsic nucleon structure, and the last terms in (4.6a) and 
(4.6b) are the inclusion of the small Darwin- Foldy correction to the radii. Note that the 
neutron (electromagnetic) contribution to the proton weak radius is not suppressed by any 
(1 — 4sin^ 9w) factor, and thus is surprisingly significant. 

Using sin^^^ ^ .23, (r^)]:^^ ^ 0.7 fm^, (r^)]^^ ^ -0.11 fm^, {r^j = gives 



r 



I,p 2 



i(l-4sin2^VK)(2.1 fm^), (4.7a) 



r^Y ^ -i(.74fm^), (4.7b) 



I,n 2 

The quantities in parentheses above can be interpreted as the physical size (squared) of the 
weak distributions. Note that using numbers of 0{±.l) |12[ for the strangeness radius will 
have a large effect on (r^)^ for both the proton and (somewhat less so) the neutron. 

To a good approximation, considering only rms radii, but no higher moments, the relevant 
PNC matrix element is then given by a convolution of (point) nucleon centers with their 
intrinsic structure, yielding a replacement for Eqs. (2.3) and (2.8), 



<i\HpNc,i\j> 
G p 

^71 



< 



(ivg>„(f ) + zg>,(f ))v^l7^^e> d'r. 



(4.8a) 



^^C,,{Z)M[NQlqr. + ZQ-^q.l (4.8b) 



with the quantities and g„ slightly modified from eqns. (2.9), 



(4.9) 



where Pp„(r) is now the density distribution of nucleon centers, normalized to 1. 
Assuming, for simplicity, uniform nucleon distributions, with Rp ^ R, 

<lp ^ l-(^«)^(.26 + 1^(2.1 -(r^);/2Q-)), 
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(4.10) 

with all radii measured in fm. 

For ^°^Pb, the internal nucleon structure contributes about 0.002 to g„, and a possible 
strangeness radius discussed above (0.14 fm^) would contribute about 5 times less. The 
internal structure corrects qp by about 0.005, and the strangeness contribution here would 
be comparable, about 0.004. In Cs, these numbers turn out to be smaller by about 40%. 

Prom the discussion to come in Sec. VI, we will see that these contributions from 
(known) finite nucleon structure contributes at about the 0.2% level in an extraction of the 
weak nuclear charge when measuring a single isotope of Pb (0.1% level for Cs). This might 
need to be taken into account in an extremely high precision analysis, but it will not add to 
the uncertainty in testing the Standard Model. (See also the complete discussion in Sec. VI 
to compare with the expected scale of nuclear, atomic, and electroweak radiative corrections 
and uncertainties.) On the other hand, strangeness contributions, which are currently very 
uncertain, might affect a determination of the weak charge at below the 0.1% level in Pb, 
and even less in Cs, and thus are likely to be quite negligible. They could only become 
relevant if the nucleon strangeness radius were comparable to the electromagnetic radius 
itself. 

In the case of isotopic ratios, the internal nucleon structure plays an even smaller role. 
This is because errors then come from uncertainties in the difference — Qn (see Sec. VI). 
To a good approximation, nucleon structure effects are simply additive in mean square radii 
and thus cancel in the differences. Thus neither nucleon structure, nor the uncertainties 
therein, are significant when extracting sin^ 9w from isotope ratios. 
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V. NUCLEAR MODELING 



From the rather simphstic model of section III, we aheady observe that the desired high 
precision measurements of electroweak parameters will require knowledge of neutron radii 
in heavy nuclei to within at least several percent, (see also the discussion in section VI) 
At this level, one clearly must treat higher moments with some care, and the microscopic 
details of the nucleon distributions may be of some importance. For this reason, we have 
attempted to evaluate and Qp numerically, utilizing the best existing nuclear models for 
neutron, proton, and charge distributions available to us. In this section, we discuss some of 
the basic features of these models, along with some caveats on their reliability for neutron 
observables. 

The nuclear many-body problem presents a formidable challenge for infinite nuclear 
matter, and an even greater one for heavy finite nuclei. The most popular route being taken 
today is some version of Hartree-Fock, which has had considerable success in describing a 
variety of nuclear properties semiquantitatively. 



A. Brueckner-Hartree-Fock 

The underlying basis of nuclear Hartree-Fock calculations is Brueckner theory. The 
elementary two-body interactions are too strong (especially the short-range repulsion) to 
lead to meaningful HF calculations. Although there has been extensive work on nuclear 
matter calculations using Brueckner theory and beyond, for finite nuclei only light ones 
have been considered |13| and nothing for the nuclei of interest here. 



The lowest level is the independent pair approximation. The effective interaction is not 
V but the Brueckner G-matrix, where G = vF, with ?/'(!, 2) = F(l, 2)0(1)0(2). G satisfies 
a scattering-type equation with a projection operator in intermediate states which excludes 
scattering back into the Fermi sea; this also results in no phase shift (for pairs in occupied 
states) and to the two-body wave function "wound" which extends over a "healing distance." 
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The 0's are to be identified with the HF single guasi-particle functions. We have used F 
here to denote the two-particle correlation function. For repulsive core potentials, F has a 
hole (wound) centered about r = 0. The Brueckner G is non-local, and both energy- and 
density-dependent . 



B. Two-body correlations 



The dependence of the single particle density distribution on the correlation function is 
relatively small. We can estimate it as follows. Let 



P2(fi, fa) = Pi(ri)pi(r2)/(fi - fa) . 



(5.1) 



Let pi(r) oc e '"^/"^ and /(r) oc e For << oat, one finds for the rms size of the single 
particle density distribution 



j d^ri d^r2 rl p(fi,f2) 



1 ((hV 
4^2 Vajv/ 



(5.2) 



where <r^>i corresponds to pi. For heavy nuclei (say aN/dc ~ 7.0/0.7 fm), the correction 
is less than 2 x lO"'^, which is below our level of current concern. 



C. Phenomenological Hartree-Fock, including deformation and pairing 

Because of the numerical complexity, most HF calculations have employed phenomeno- 
logical potentials intended to simulate the Brueckner G-matrix. The most commonly used 
potentials are varieties of the very convenient Skyrme interaction. The Skyrme interactions 
are of the delta-function form and as such lead to single particle equations with local one- 
body potentials and spatially-dependent effective masses, with no more complication than 
Hartree calculations. In contrast, finite range interactions lead to non-local single particle 
potentials arising from the exchange term. Momentum-dependent Skyrme interactions do 
not lead to further complications and simulate some effects of finite range. Calculations 
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have also been done with finite range forces, using e.g. the Gogny interaction |TJ] Note that 
none of these phenomenological potentials are intended to reproduce free nucleon-nucleon 
scattering. There are of the order of eight (more or less) adjustable parameters in any model. 



Most nuclear structure calculations on heavy nuclei are carried out in the deformed 
Hartree-Fock or the Hartree-Fock-Bogolyubov approximations. The latter include BCS-type 
pairing. Hartree-Fock encompasses a limited class of correlation structure. Only correlations 
of a collective nature are included. It is not surprising to find that in HF neutron and proton 
densities tend to track one another. Nevertheless, they do exhibit the expected behavior 
that the neutron rms radius increases more rapidly than the proton one with increasing A 
in an isotopic sequence. The relative tracking (variation in the neutron skin) depends on 
the way in which symmetry energy is handled. 

Intrinsic deformations play a key role in spherically averaged proton and neutron densi- 
ties. In the uniform, incompressible approximation, for example, the mean square radius is 
increased according to PB[ 



1 + ^<P'> 

4:71 



(5.3) 



where j3 is the nuclear shape parameter, proportional to the quadrupole moment. (3 can 
attain values of the order of 1/3 and changes in < > among isotopes can produce deviations 
in spectroscopic isotope shifts by an order of magnitude from the A^^^ law. Although HF 
calculations tend to yield spherical {<(3 = 0>) or near spherical equilibrium shapes for the 
Pb isotopes, < > is not zero and changes in < > have been considered by some authors. 



D. Relativistic Mean Field 



Although nuclear structure is primarily nonrelativistic, considerable success has been 
achieved by treating the nucleons and protons as point Dirac particles [|T7| [see, however, 
Achtzehnter and Wilets [|TB[] interacting with phenomenological vector and scalar mesons in 
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the mean field approximation. The vector mesons can be identified with the isoscalar omega 
and the isovector rho mesons; the scalar meson is a simulation of two-pion exchange. The 
mesons are treated in the mean field, or c-number approximation. An attractive feature of 
the model is that the strong spin-orbit potential appears to emerge "naturally." While this 
turns out to be true for isoscalar potentials, it fails badly for the isovector potentials, but 
can always be parameterized to yield reasonable results, [l^ 

In order to fit nuclear properties, it has been necessary to go beyond linear field theory. 
Self- interact ion of the scalar field has been introduced, with additional parameters. Among 
other problems, this solved the compression modulus anomaly, which is much too large in 
the linear model. The total number of adjustable parameters which must be introduced is 
comparable to that required in models using Skyrme forces. As with Skyrme forces, the 
mean field approximation does not lead to nonlocality in the one-body potentials. 



E. RPA and MCHF 

The tail of the neutron or proton distribution has correlation/polarization corrections 
not described by HF, at least for large distances. The reason for this is that the individual 
nucleon wave functions see the potential of the "mean" self-consistent core. In the tail 
region, the residual core tends to relax. This is most evident for the separation energy: In 
HF, the separation energy of a nucleon is just the energy eigenvalue (Koopmanns' theorem) 
if the core is frozen. If the energy of the residual nucleus is recalculated self consistently, the 
separation energy is reduced by what is termed the rearrangement energy. Correlations of 
this type are included through RPA, which is equivalent to small amplitude, time- dependent 
Hartree-Fock. 

Other types of correlations could be included through multiconfigurational Hartree-Fock, 
which, as the name implies, means that the trial wave function is not a single determinant, 
but a sum of determinants (configurations). This serves two purposes: correlations of the 
kind allowed by the choice of configurations are included, and the occupation of these con- 
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figurations modifies tlie mean field potential and single particle functions. 



F. Beyond Hartree-Fock plus 

Most HF and HFB calculations do reasonably well in reproducing atomic isotopic shifts 
for the even-even isotopes of PB below 208. So, incidentally, does the droplet model of 
Meyers [|l^]. They all fail to reproduce even-odd staggering, which shows odd nuclei to be 



smaller than the mean of their even-even neighbors, and also do badly on the shifts above 
208. 

There are no giant shell model diagonalization calculations available which yield densities 
for heavy nuclei. Such would be very valuable for comparison with Hartree-Fock results, 
since in principle they include all types of correlations, limited only by the size of the basis. 

An idealized shell model calculation should be based on realistic two body interactions, 
the kind which fit free two-body scattering data. The Hilbert space could be divided into a 
"near" and a "far" space. The effective two body interaction could be obtained by solving 
for the Brueckner G-matrix with the intermediate states excluded from the near space. The 
far space scattering states could be approximated by plane waves if the momentum sphere 
separating the spaces is sufficiently large. |]l3l The Hamiltonian matrix for the inner space, 
using the effective interaction, is then diagonalized. 



G. Summary and Discussion 

Unfortunately, not all of the theoretical considerations discussed above have been incor- 
porated in any single calculation. Heavy nuclei pose a difficult challenge for reliable, detailed 
modeling at the level of precision we require. There do exist in the literature a number of 
recent efforts, as discussed in sections C and D above, which involve either relativistic or non 
relativistic Hartree-Fock nuclear calculations. We have accumulated densities from several 
of these authors in order to evaluate qp and g„ and make comparisons among the different 
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models. These include various HF calculations with Skyrme forces, pO|,pT| an HFB calcu- 



lation with a Gogny finite range DIS interaction, |I^,p2| and several relativistic mean field 



models. p3| , p4| The results are presented in Sec. VII. 

The modelers fit their adjustable parameters to choices among various bulk properties 
(energy per nucleon, compressibility modulus, symmetry energy, etc.), and properties of 
particular nuclei (energies, charge radii, deformations, spectra, multipole sum rules, etc.) 
Indeed, the physics behind the models comes, in part, from the choice of the particular 
observables included in the parameter fits. The models we have selected all do roughly 
equally well in fitting the wide range of nuclear observables available across the periodic 
table. [|l],|3l 

The analysis of atomic PNC, as discussed in sections II and III, relies on a detailed 
knowledge of neutron distributions in nuclei. The lack of unambiguous, precise experimental 
measures of neutron radii means that all of these models must "extrapolate" to the desired 
neutron properties. Charge radii, on the other hand, are in a certain sense "built in" , in that 
the set of observables to which the nuclear model parameters are fit includes charge radii 
of several even-even nuclei, one of which is ^''^Pb. In defense of the models, they predict 
with good success the charge radii of other even-even nuclei not included in the fit, and 



also reproduce the well-measured isotopic charge radius shifts for, e.g. the even isotopes of 
Pb. p0|j37| However, they do not reproduce the observed even-odd staggering of charge radii 
very well, nor are the results as good for the charge radius of ^^°Pb, an indication that some 
care should be taken when considering non closed-shell cases. 

There do exist some data which may give more direct information on the neutron skin. 
This might be used as additional input to these nuclear models, and could help further con- 
strain the predictions for neutron radius, and neutron isotope shifts, if one could demonstrate 
consistency in the results. Perhaps the best known data comes from 800 MeV polarized pro- 
ton scattering from ^o^Pb. [^,^ This gives Rn — -Rp = 0.14 ± 0.04 fm. The quoted error, 
which is quite small for our purposes, contains both statistical and certain theoretical uncer- 
tainties as stated. However, there are still additional theoretical uncertainties, involving e.g. 
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assumptions about the in medium nucleon-nucleon t-matrix, and the result seems to exhibit 



a rather large and troubling energy dependence. The absolute value of i?„ in one isotope 
is believed to be fairly difficult to obtain with confidence using such experiments. However, 
it may be that the relative shift in Rn among isotopes involves cancellations that reduce 



these theoretical uncertainties. Measurements on other Pb isotopes (data apparently 
exists for ^°^Pb) would clearly be of interest in this context. Further experimentation and 
theoretical analysis at other energies are also crucial to demonstrating the consistency of 
the results. 

Experiments involving intermediate energy charged pion scattering from nuclei may also 
help further constrain the neutron radii, or the relevant isovector model parameters, as well. 
Such data exists for ^°^Pb, [0 but again the absolute normalization poses a real challenge 



to analyses. Taken at face value, the p and vr results for lead neutron radii do agree 



with one another reasonably well, and also match with e.g. the Gogny finite range Hartree- 
Fock calculations. Another experimental possibility involves energies and sum rule strengths 
of giant multipole resonances. The uncertainties here are even larger, and difficult to 
estimate. Clearly, a reliable set of such additional "strong probe" inputs, including yet other 
options such as a particle scattering, kaon scattering. Coulomb displacement energies, etc., 
could be an aid in constraining the theoretical models on neutron properties. 

Another promising experimental possibility for the future might be direct electroweak 
experiments, such as parity violating asymmetries in elastic, intermediate energy e~ - nucleus 



scattering, as proposed by Donnelly, Dubach, and Sick, [^] or perhaps elastic u scattering. 
The reactions and analyses are quite clean, just as in the charge scattering case. There would 
be, for example, no serious ambiguities in the absolute scale of the radii measured. Such 
experiments would in fact be sensitive to the full nuclear weak charge distribution, rather 
than just the RMS radius, which many of the strong probe measurements are primarily 
sensitive to. Because such experiments could be done at moderate momentum transfers 
(g ~ 1 fm^^), the extraction of nuclear distribution information would be much less sensitive 
to the precise values of electroweak parameters than in the corresponding atomic parity 
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violation case. The asymmetries and v cross sections are naturally extremely small, and 
the experimental challenges are formidable. Nevertheless, recent estimates for the parity 
violating asymmetries indicate that measurements sensitive e.g. to the neutron RMS radius 
in 208p]^ ^]-^g level are feasible. As we will see in sections VI and VII, such a level 
would make the nuclear structure uncertainties quite negligible for the purposes of extracting 
standard model parameters from single isotope atomic PNC measurements. 

In any case, current model fitting has been done with the best and most reliable data at 
hand, most of which are not directly sensitive to neutron distributions. It is always difficult 
to estimate the theoretical uncertainties in such model calculations. In this section, we have 
already mentioned several potentially important missing features that future work should 
address, especially involving nucleon correlations. We have not attempted here to try to 
choose a "best model" from the various ones we examined, but rather wish to evaluate the 
existing spread in predictions as an effective lower bound on the theoretical uncertainties 
involved. One might, however, try to make a selection based on detailed comparisons, specif- 
ically targeting a good fit to heavier nuclei energies, isotopic shifts, giant dipole properties, 
and other quantities potentially sensitive to isovector properties. We encourage work in such 
directions. The goal should be to find the most reliable model(s) while still retaining an 
estimate of the remaining theoretical uncertainties. 



VI. ERROR ANALYSIS AND TESTS OF ELECTROWEAK PHYSICS 

One of the motivations for further improving atomic parity violation experiments is 
to test the Standard Model at the level of its one-loop electroweak radiative corrections. 
This allows one to probe for possible small "new physics" effects, which would appear as 
further loop corrections or more directly as additional interactions at the tree level (i. e. 
without loop corrections). A good example of the latter, to which atomic PNC is particularly 
sensitive, is exchange of a second, more massive, neutral Z-boson required in theories with 
larger gauge groups. In the Standard Model, the loop contributions are separated into two 
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parts: fixed radiative corrections due to contributions from tlie known quarks, leptons and 
bosons, and tlie lieavy pliysics part due to contributions from the top quark and the Higgs 
boson. One is interested in an experimental determination of the heavy physics part, which 
in the language of Marciano and Rosner |^ is expressed in terms of weak isospin-conserving, 
S, and isospin-breaking, T, effects. These two constants ||2^ are a convenient way not only 
of including uncertainties in the top quark and Higgs masses, but of parameterizing the 
effects of some specific classes of new physics as well. It turns out that low energy PNC 
measurements are nicely complementary to high energy measurements such as direct Z- 
boson production, since both the radiative corrections and the sensitivities to new tree-level 
interactions are quite different. 

To show how the radiative corrections and the possible new physics enter into atomic 
PNC, and how they might compare in size to nuclear structure effects, we begin by rewriting 
equation (2.9) in the form: 

<t\Hp^c,i\j>= ^Q,iZ)Af[QwiN,Z) + QZ''{N,Z)] (6.1) 

where Qw{N, Z), known as the nuclear weak charge, is the quantity of primary interest to 
electroweak theory, and in the standard model without radiative corrections reduces to: 

Ql, = -N + Z{l-4:x) (6.2) 

where x = sin^ 9w- Qw{N, Z) is determined from atomic experiments by combining atomic 
measurements of < i\HpNc,i\j > with calculations of both atomic structure (contained in 
the factor Cij) and nuclear structure. The nuclear structure corrections are contained in 
Qw^i^' Z), which is given by: 

Q'^'{N,Z) = Qw{[qn-l]N,[qp-l]Z) 

^ -iV(g„-l) + Z(l-4x)(gp-l). (6.3) 

Nuclear structure is also contained in the normalization A/", but as we will see in section VII, 
A/" is determined by the nuclear charge distribution, which is usually known experimentally. 
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When we include possible new physics, together with the effects of radiative corrections 
which have been calculated by others, ||2^ Qw{,N^ Z) becomes: 



Qw{N, Z) = (0.9857 ± 0.0004)(1 + 0.00782T) 
X [-A^ + Z{1 - (4.012 ± 0.010)^] 

+Ql7e{N,Z) (6.4) 

where x is assumed here to be defined at the mass scale mz by modified minimal subtraction. 



25| and is given by: 

X = .2323 ± .0007 + .003655 - .00261T . (6.5) 

The errors indicated in (6.4) and (6.5) come from uncertainties both in experimental input 
parameters and in evaluations of known physics loop-diagrams. The unknown, heavy physics 
loop-corrections are contained in the parameters S and T, which depend upon the heavy 
masses, and are defined such that 5 = T = if rriH = 100 GeV, rrit = 140 GeV, and if there 
is no new physics beyond the Standard Model. Including Q^^^{N, Z) in Qw{N,Z) allows 
for additional tree-level physics beyond the Standard Model. For example, exchange of the 
extra Z^ in SO(IO) models p5| , p7i (assuming no Z^ — Z mixing) would make: 

QZ^,{N, Z) ^ 0.4(2iV + Z)ml,/m\ . (6.6) 

It is useful to consider how well the parameters in Qw are currently known. The central 
value of Q^^^{N, Z), determined mainly by Cs PNC measurements, is about 2.2 ± 1.6 ± .9 
{if all other heavy physics in equation (6.4) is ignored), and corresponds in the SO(IO) 
model to niz^ ~ 500 GeV. Conversely, assuming no new tree-level physics (i.e., Q^/Je — 0)i 
the experimental uncertainty in T is currently around ±1, and in S around ±3, the latter 
determined largely from Cs PNC. Ultimately, as Marciano and Rosner have indicated, an 
effort to reduce the uncertainty in S to ±0.2 is extremely important, since at that level it is 
sensitive even to minimal one-doublet technicolor models. This sort of accuracy is an extreme 
challenge to either high energy or atomic experiments. Current knowledge of sin^ 9\y from 
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a global analysis of electroweak data |2^can be summarized hj x = 0.230 ± .004, (roughly 



2% uncertainty). If future high energy measurements were to reduce the uncertainty in 
X beyond what is attainable in atoms, the atomic experiments would still be valuable for 
improving the limits, e.g., on an additional Z. 

In summary, any improvement in determining atomic PNC is likely to provide useful 
information about electroweak physics, and it becomes extremely important to work out 
how much nuclear structure uncertainties may be a limiting factor, and to reduce these 
uncertainties where possible. 

We first consider the impact of nuclear uncertainties on PNC measurements of single 
isotopes. PNC experiments to date have been done on stable isotopes of heavy atoms, namely 
Cs, Pb, Bi, and Tl, and have not compared different isotopes of the same element. From 
equation (6.1) we derive an expression for the uncertainty in Qvy in terms of the uncertainties 
in atomic and nuclear structure and in the measured quantity O =<i\HpNc,i\j >'■ 



nuc 



6Qw SO 6Cij 6N' 



(6.7) 



Qw ^ M Qw 

If we assume that O can be measured to arbitrary accuracy, and that proton distributions 
(which will influence dM) are also well enough understood and/or measured, there remain 
the uncertainties coming from atomic and nuclear structure, which we can write in the form: 



5Qw SCij 



- (6.8) 



Qw Cij 

where we have dropped all terms containing the factor 1 — 4.0 12x, which should be quite 
negligible due to the accidental value of x ~ |. Rewriting in terms of the weak interaction 
parameters, we obtain: 

iV C^iy L^ij 

ignoring the contribution of the weak-isospin breaking parameter T which cancels to better 
than 10% for the full range of Z/N found in the elements of experimental interest. Thus a 
PNC measurement in a single isotope can set limits on the weak-isospin conserving parameter 

26 



S and/or new tree- level interactions, and in fact the best limits on both of these parameters 
now come from PNC measurements in atomic cesium. To determine the role of nuclear 
structure, we must compare the uncertainty 6qn on the right hand side of equation (6.9) to 
the atomic structure uncertainty 5Cij/Cij . This we do later, in section VII. 

Because of the difficult atomic physics calculations, there has been some serious interest 
in measuring parity violation in a chain of isotopes. Taking ratios between isotopes cancels 
essentially all dependence on atomic structure. Unfortunately, although the atomic physics 
indeed cancels in the ratio, the nuclear structure does not. Referring to equation (6.1) we 
consider the ratio: 

_ [QwjN, z) + (iV, Z)]H 

O' [Qw{N', Z) + Q^'^iN', Z)\M' ■ ^ ■ ^ 

where primed and unprimed quantities refer to different isotopes. The sensitivity of x and 
Qtreei-^1 extracted from this ratio, to the nuclear structure is then given approximately 

by 



new 



NN'lSn S(AM) 5fAoJl , ^, 

(6.11) 



sn ^ 5iA^^) ^ siAq^) 



UN Qn 

where we have made simplifying assumptions that the isotopes are close together, i.e. AA^ = 
N' — N « N, that sin^^vi^ fa 1/4, and where we have used e.g. {5qn/qn — ^Qn/Qn) ~ 
5{Aqn)/qn, which is numerically accurate for the models of Pb we have considered. Because 
of the special sensitivity of atomic PNC to any additional heavy Z-bosons, we note as an 
example that a determination of rriz^ in the model of equation (6.6) would be constrained 
by replacing the left hand side of equation (6.11) by [zbxjx — 0AZ5(m^^/m%J^. 

The uncertainties on the right side of (6.11) are effectively in the relative difference 
between quantities for two isotopes. In principle, different nuclear models which disagree on 
the absolute values of, say, qn may agree on the relative change in this quantity to a much 
higher degree of accuracy. However, such a reduction in uncertainty in the terms within the 
brackets in expression (6.11) is roughly compensated by the factor N/AA. Comparing with 
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equation (6.9) for a single isotope, in which any new tree-level interactions enter with the 
equally uncertain loop-parameter S, we see that when we instead compare isotopes, Q^/Je 
appears together with a different parameter, x, which is independently measurable in high 
energy experiments. 



VII. DISCUSSION 

To calculate we use various theoretical predictions of neutron and proton distribu- 
tions from the literature. Proton distributions are used to compute f{r), the electronic wave 
function overlap defined in equation (2.5). We solve numerically for single electron Dirac 
si/2 and pi/2 wave functions near the origin, in the Coulomb potential of the nuclear charge 
distribution (as discussed in section II), and make no approximation of a power series in 
Za, as was done e.g. for equation (3.2). We have neglected the contributions to the nuclear 
charge distribution from internal neutron structure, as discussed in section IV. We estimate 
the error associated with this assumption to be well below the level of the model uncer- 
tainties themselves. The quantity Qn = J Pn{r)f{r)cPr is then calculated directly from the 
corresponding neutron distribution. 

In Tables | and H, we present the rms radii Rn and Rp, the correction factors g„ and g^, 
and the electron normalization A/", for several nuclear models of the Pb isotopes 202 and 210. 
Except for the norm, the spread in values in the final rows of Tables | and |T| should give some 
indication of a lower limit on the current level of theoretical model-dependent uncertainties, 
assuming that one accepts these models as equally phenomenologically reasonable. 

The normalization factor Af defined in equation (2.5), which is proportional to l/f{r ^ 
Rp), is defined arbitrarily here as 0.10361//(300 fm). The numerator, /''^p*(300 fm), is eval- 
uated using a model independent experimental charge distribution from electron scattering 
off ^°^Pb. As stated earlier, we are not concerned with the absolute value of the norm, but 
only its dependence on atomic weight and charge distribution. This definition simply scales 
A/" to be near 1.0. For ^°^Pb, the model spread in the normalization from Table ^ might 
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appear to contribute at a significant level. One can, however, consider correcting TV by using 
an approximate formula relating M to the charge radius, namely M = R~'^ , as in eqn. (2.6). 
This is given by 

m'=u{rvr:^^,). (7.1) 

The model spread in this M' is significantly reduced. The point is that these models are 
not precisely reproducing the observed charge radii of the lead isotopes, which feeds rather 
directly into a calculation of M . The correction factor above compensates for this, using the 
existing high precision measurements of charge radii from optical isotope shifts and electron 



scattering. 37,3 



Some of the results in tables | and || are reproduced in graphical form also, in Figures 
m, H, and 1^. In Fig. |l], we plot the predicted versus atomic weight for several even 
lead isotopes, and in Fig. ^ the ratio Rn/Rp. The spread in i?„ among models is decidedly 
larger than the spread in Rp. We do note a systematically larger neutron radius in the 
relativistic models. The origin of this is indeed not yet completely understood, but may 



be connected with larger asymmetry energies found in these models. This in turn might 
tend to pull neutron and proton distributions together where the densities are high, leaving 
a somewhat larger neutron tail. 

In Fig. ^, we plot g„ versus atomic weight for several even lead isotopes. The spread 
is closely related to the spread in Rn/Rp shown in Fig. |^, as might be expected from the 
simplified formulas (1.1) or (3.3b) based on uniform nuclear charge density. Estimates of g„ 
using these simplified formulas yield the same general trends as the detailed calculations, 
with absolute values differing generally by parts in a thousand or less. The relativistic 
models yield somewhat smaller due to their larger R^/Rp ratio. 

In the case of single isotopes, the total nuclear model spread does not appear to be the 
most serious problem in using equation (6.8) or (6.9) to extract weak interaction parameters 
from atomic PNC . For ^°'^Pb, the typical full spread in calculated g„ is ^ 0.005. In the case 
of ^^'^Cs, the sensitivity to nuclear structure is even weaker, due to the smaller value of (Za). 
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A larger uncertainty, at least at the present time, is due to atomic physics calculations. 
For example, Cs is one of the most favorable elements from the point of view of atomic 
theory, and to achieve the current level of quoted uncertainty of 5Cij{Z) /Cij{Z) ^1% in 
Cs is an impressive task. But this uncertainty is still probably larger than the uncertainty 
in g„ for Cs. Significant future improvement in atomic calculations is likely to be difficult. 
Thus, aside from any experimental uncertainties, atomic structure is the present limiting 
factor in getting Qw and the associated weak parameters from single isotope atomic PNC 
measurements, and appears to remain so even after considering the possible nuclear physics 
effects. This conclusion is consistent with the findings mentioned in the calculation of 
reference 2. 

Consider next the ratios in an isotopic chain, for example (^°^Pb/^°^Pb). Referring 
to equation (6.11) we see that the PNC experiments would then be measuring x and/or 
observing new tree level physics. For definiteness let us assume no new tree level physics. 
Then a ±1% extraction of x would require (5Ag„ ^ 6 ■ lO"'^, and 6AAf/Af ^ 6 ■ 10~^. 
Assuming uniform nuclear distributions, this implies 6A{Rn/Rp) ^ 4 ■ 10~^. Referring to 
Table |T|, which shows the change in various quantities between these two particular lead 
isotopes, the model spread for (5Ag„ is around 9 ■ 10^"^ and for 6AAf/Af is about 6 ■ 10~^. 
Note however that when JV is corrected as in (7.1) above, using experimental knowledge of 
charge radii, this spread, at least, is significantly reduced to below the ibl% level. This is 
seen from the final column in Table |T|. But the spread in Ag„ remains, and is comparable 
to the accuracy needed for a 1% extraction of x. Similarly, the model spread in (Rn/Rp) 
from Table III is about 5-10"'^, which likewise corresponds to a >1% spread in x. 

If we exclude the relativistic models, which seem to have substantially different neutron 
radii from the conventional H-F calculations, the model spread just among the various 
Skyrme parameterizations considered gives (5Ag„ ^ 4 ■ 10~^. It thus appears unlikely that 
PNC measurements comparing Pb isotopes could yield much better than a 1% determination 
of X, unless there is significant improvement in understanding of nuclear structure. 

The same results can be seen perhaps more clearly in Fig. which displays in graphical 
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form the values of Ag„(202 208) from Table |T| versus the different models considered. 
The spread in predictions of this quantity is actually larger than 100%. Also shown in the 
figure is a typical scale of 1% in the weak angle. As noted already, the model spread is too 
large for extractions of x at the sub 1% level if one cannot otherwise eliminate or improve 
any of the models used. On the other hand, the nuclear structure uncertainties may not 
preclude a significant improvement in sensitivity to new Z bosons or other new tree-level 
physics in equations (6.11), particularly if x is determined well by high energy experiments. 

Although the nonrelativistic models do appear to cluster together somewhat, one should 
perhaps be a bit wary of their apparent self-consistency. For example, a modification of 
the coefficient of the isovector (n-p asymmetry) surface term, a (pp — p„)V^(pp — p„) term 
in the Skyrme Lagrangian, ||21| has little effect on most bulk properties, and hence on 



the goodness of the Skyrme fits. This term, however, does modify the neutron skin 
significantly. Reinhard's rough estimates show that an uncertainty of Ai?„ ^ ±.15 fm is 
not unreasonable. This in turn can modify the quantity shown in Fig. ^ by amounts of 
0{6 ■ 10~^), larger than the spread in the given Skyrme models. The relativistic models do 
not have such fiexibility, as the isovector rho couplings are largely constrained by isotopic 
trends in ground state energies and charge radii, but this is of course no guarantee that 
these models correctly describe all isovector properties equally well. 

In the case of Cs isotopes, accurate calculations for neutron radii (or even proton radii) 
are difficult. They have odd Z, and require additional approximations to deal with unfilled 
shells, as well as deformations. The lack of success in predicting the even-odd staggering of 
5{f'^)ch in lead isotopes indicates the seriousness of these problems. An estimate of the scales 
involved, however, can be made using calculations with existing nuclear codes. One such 
result g gives Rn/ Rp ~ 1.03 for ^^^Cs, and AR^jRp (^^^Cs ^ ^^^Cs) ^ 5 ■ 10~^. If this 
latter number itself has a 100% uncertainty (for comparison, see Fig. ^ for the case of lead 
which does show a 100% spread among model predictions of the equivalent quantity Ag„ for 
about the same AA/A), then the uncertainty in x from this fairly small range of isotopes 
would be approximately 1%. AA of up to 10 or higher may be experimentally possible for 
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Cs, which might help to reduce the nuclear physics uncertainties. Prom the experimental 
side, the absence of stable isotopic partners to ^^^Cs makes it difficult to obtain values of 
ARn/ Rp from parity violating electron scattering, or p elastic, or pion experiments, as may 
be possible for the lead isotopes. Further work on theoretical estimates for Cs isotopic radii 
is clearly called for. 

Given a set of experimental results for isotopic PNC ratios, one can also consider a boot- 
strap procedure: from atomic experiments over several isotope differences, use the various 
models to extract the weak mixing angle. Then, only those models which yield the same 
sin^ 9w for the various isotopic pairs are acceptable. Unfortunately, the various nuclear mod- 
els we have considered (for lead isotopes near ^°^P6) yield predictions for the PNC ratios 
which are fairly linear with A^l. Since this prediction is also roughly linear with sin^ 9^, 
it appears that the various nuclear models could be internally consistent, each yielding a 
unique sin^ 9w but differing from model to model about the extracted value. Of course, 
one cannot draw any firm conclusions about this until after the data are known. There are 
indeed some slight deviations from linearity, especially for non-closed shell isotopes, and one 
may be able to take advantage of this. In essence, this bootstrap idea uses PNC atomic 
isotope ratios themselves as our desired additional constraint on neutron properties - with 
a large enough set of PNC data, one could hope to simultaneously constrain the nuclear 
model parameters and measure the weak mixing angle. 

VIII. CONCLUSIONS 

For the case of Pb, in order to extract electroweak parameters from atomic PNC ex- 
periments at a level of precision which would be considered "significant" for testing the 
Standard Model, we have shown that it is necessary to have confidence in the isotopic rel- 
ative neutron/proton radius shift, A{Rn/ Rp)/ AA, to better than a few times 10""^. We 
have examined various nuclear model calculations, and find that the spread in theoretical 
values corresponds to an uncertainty in the weak mixing angle greater than 1%, with the 
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assumption that no new physics is present. Without some further basis for discriminating 
among the various models, the spread represents a lower bound to the uncertainties in the 
calculated values. 

The basic problem is essentially that the models have been parameterized to fit proper- 
ties like charge distributions, which are not directly sensitive to neutron distributions. As 
Reinhard has shown, it appears that a surface symmetry energy term in certain nonrelativis- 
tic (Skyrme interaction) nuclear models can be "dialed" somewhat to change the neutron 
size without significantly spoiling the basic fits. Including data which are more sensitive to 
neutron properties, such as isotopic trends in ground state properties, and perhaps giant 
resonance energies and sum rules, could be useful to constrain such terms. 

There do exist experiments which are sensitive to neutron radii, e.g. tt^/tt^ scattering, 
and medium energy polarized proton scattering. If the quoted errors on the latter can be 
taken literally, one could use it to discriminate among the various models and provide the 
confidence one needs to extract the desired electroweak parameters from atomic experiments. 
It would be valuable to repeat the experiments and analyses at other energies in order to 
demonstrate the consistency of the results, and to consider both vr^ and p scattering on 
multiple Pb isotopes for a direct experimental measure of the isotopic shift in neutron radii. 
We have also noted in this work that the detailed distribution of neutrons, beyond just the 
RMS radius, is of some importance. This implies that we may still have to rely on the 
nuclear models for an extraction of the electroweak parameters. As discussed earlier, the 
use of alternative electroweak probes, such as parity violating (polarized) electron scattering 
at intermediate energies, |^ would be of obvious value for independently extracting the 



desired neutron distribution. 

We can turn the problem around, however, and note that an accurate measurement 
of X from high energy experiments presents a unique opportunity to extract the isotopic 
neutron radius shifts from atomic experiments cleanly, and hence test the nuclear models. 
The situation is quite analogous to the extraction of changes in charge radii from atomic 
isotope shifts. 
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FIGURES 

FIG. 1. Root mean square neutron radius, in fermis, for lead isotopes, plotted versus atomic 
weight. Models are defined as in Table |. We also show points for several additional relativistic 



Hartree parameterizations, NLl, NL06, and NL075. |23,24| Skyrme calculations are connected with 
dashes to guide the eye. Gogny HFB is connected with dots. Relativistic models are connected 
with dot-dashes. 

FIG. 2. Ratio of neutron to proton radius for lead isotopes, plotted versus atomic weight. 
Symbols are defined as in Table | and Fig. 

FIG. 3. Neutron correction factor g„ for lead isotopes, plotted versus atomic weight. Symbols 
are defined in Table | and Fig. |l|. The large error bar on the right side represents the allowed 
spread corresponding to a uncertainty in the weak charge, as given by eqn. (6.8). The 

vertical position of this error bar is arbitrary. 

FIG. 4. Change in Qn between '^'^^Pb and ■^'^^Pb, shown versus model weight. (The x axis is 
arbitrary) The error bar on the right side represents the allowed spread in the plotted quantity 
corresponding to a ±1% uncertainty in x, as discussed following eqn. (6.11). The vertical position 
of this error bar is again arbitrary. 
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TABLES 

TABLE L Some properties of ^''■^Pb relevant to atomic parity violation, for several nuclear 
models. Properties listed are r.m.s proton radius, neutron radius, ratio (Rn/Rp), difference in 
mean-square charge radius from ^'^^Pb, qp and g„, defined in eqns. (2.9a and b), the normalization 
factor, M, and a "renormalized norm" defined in (7.1). All distances are in fm. The models listed 
are Hartree-Fock-Bogolyubov with a Gogny ||2^ finite-ranged DIS interaction, using a 15 shell 
spherical harmonic oscillator basis (G:HFB), various parameterizations of the Skyrme interaction 
pO| , ^ in the spherical Hartree-Fock approximation, Skyrme A (SkA), star (Sk*), and 3 (Sk3), and 



relativistic Hartree mean field calculations : 24 ] with a nonlinear PL40 parameter set (Rel) . The 



first row contains experimental numbers, where known. [37,38,4^] The final row simply indicates 
the maximum spread among the models. 



202pb 






Rn/Rp 




qp 




M 


M' 


Data 








-.330(4) 










G:HFB 


5.409 


5.519 


1.020 


-.274 


.90599 


.90318 


1.0087 


1.0041 


SkA 


5.431 


5.607 


1.032 


-.297 


.90604 


.90141 


1.0069 


1.0039 


Sk* 


5.423 


5.560 


1.025 


-.296 


.90608 


.90254 


1.0076 


1.0039 


Sk3 


5.488 


5.590 


1.019 


-.355 


.90638 


.90382 


1.0018 


1.0028 


Rel 


5.484 


5.742 


1.047 


-.319 


.90618 


.89949 


1.0023 


1.0034 


Spread~ 


.08 


.2 


.03 


.08 


.0004 


.004 


0.007 


.001 
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TABLE II. Same as Table I, for ^ospb. 



208 TDK 


VP 








Qv 


(In 


Al 


Al 


Data 


5.453(/j 


5.59(4j 


l.L)3 


0.0 


.906(1) 




=1.0 




G:HFB 


5.435 


5.569 


1.025 


0.00 


.90596 


.90260 


1.0068 


1.0018 


SkA 


5.459 


5.670 


1.039 


0.00 


.90599 


.90051 


1.0049 


1.0016 


Sk* 


5.451 


5.620 


1.031 


0.00 


.90605 


.90176 


1.0053 


1.0016 


Sk3 


5.521 


5.646 


1.023 


0.00 


.90636 


.90334 


0.9992 


1.0005 


Rel 


5.513 


5.822 


1.056 


0.00 


.90607 


.89813 


1.0003 


1.0012 


Spread'^ 


.1 


.3 


.03 




.0004 


.005 


0.008 


.001 



TABLE III. Changes in various nuclear properties between ° Pb and °®Pb. The models are 
the same as in Table 1. Properties listed are the change in relative neutron to proton radii, the 

change in proton and neutron correction factors g„ and qp, and the relative change in the norm. 
The final column is for the "renormalized norm" defined in eqn. (7.1). The last row shows the 
maximum spread among the models. 



202 p]^ _^ 208 p]^ 


ARn/Rp 


Aqp 




AAf 


AAf' 


G:HFB 


.0043 


3-10-5 


5.8-10-^ 


.0020 


.00227 


SkA 


.0064 


5-10-5 


9.0-10-"^ 


.0021 


.00232 


Sk* 


.0056 


3-10-5 


7.8-10-"^ 


.0022 


.00225 


Sk3 


.0040 


2-10-5 


4.8-10-^ 


.0026 


.00230 


Rel 


.0091 


i-io-'' 


1.4-10-3 


.0021 


.00218 


Spread^ 


.0051 


8-10-5 


9-10-"^ 


6-10-"^ 


1-10-"^ 
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